STATUS OF SM CALCULATIONS OFB^S TRANSITIONS 
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We report recent progress in SM calculations of b — » s transitions. We discuss the first NNLL 
prediction of the B — > X s 7 branching ratio, including important additional subtleties due to non- 
perturbative corrections and logarithmically-enhanced cut effects, and also recent results on the 
inclusive mode B — > X s £ + £~. Moreover, new results on the corresponding exclusive modes are 
reviewed. Finally, we comment on the present status of the so-called B — > Ktt puzzle in hadronic 
b — » s transitions. 
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1. Introduction 

In any viable new physics model we have 
to understand the important flavour prob- 
lem, namely why flavour-changing neutral 
currents are suppressed. Rare decays and 
CP violating observables exclusively allow an 
analysis of this problem. However, if new 
physics does not show up in flavour physics 
through large deviations, as recent experi- 
mental data indicate, the focus on theoret- 
ically clean observables within the indirect 
search for new physics is mandatory. This 
also calls for more precise SM calculations in 
the first place. 

A crucial problem in the new physics 
search within flavour physics is the opti- 
mal separation of new physics effects from 
hadronic uncertainties. It is well known that 
inclusive decay modes arc dominated by par- 
tonic contributions; non-perturbative correc- 
tions are in general rather small 1 . Also ra- 
tios of exclusive decay modes such as asym- 
metries are well suited for the new-physics 
search. Here large parts of the hadronic un- 
certainties partially cancel out; for example, 
there are CP asymmetries that are governed 
by one weak phase only; thus the hadronic 
matrix elements cancel out completely. 

Data from K and Bd physics show that 
new sources of flavour violation in s — > d 



and b — ► d are strongly constrained, while 
the possibility of sizable new contributions to 
b — > s remains open 2,3 . We also have hints 
from model building: flavour models are not 
very effective in constraining the b — > s sec- 
tor 4 . Moreover, in SUSY-GUTs the large 
mixing angle in the neutrino sector relates 
to large mixing in the right-handed b-s sec- 
tor 5 < 6 ' 7 . 

In the following we discuss recent 
progress on several b — > s observables. 



2. B 



Among the flavour-changing current pro- 
cesses, the inclusive b — > sj mode is 
still the most prominent. The stringent 
bounds obtained from this mode on various 
non-standard scenarios are a clear example 
of the importance of clean FCNC observ- 
ables in discriminating new-physics models. 
Its branching ratio has already been mea- 
sured by several independent experiments us- 
ing semi-inclusive or fully inclusive meth- 
ods s. 9 . 10 . 11 . 12 . The world average of those 
five measurements (performed by the Heavy 
Flavour Averaging Group (HFAG) 13 ) for a 
photon energy cut E 1 > 1.6 GeV reads 

B(B -» A» cxp = 

= (3.55 ± 0.24 t°;5o ± 0.03) x 10~ 4 (1) 
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Fig. 1. New NNLL prediction versus HFAG average. 



where the errors are combined statistical and 
systematic, systematic due to the extrapola- 
tion, and due to the b — > d 7 fraction. 

On the theory side, perturbative QCD 
contributions to the decay rate are dom- 
inant and lead to large logarithms a s x 
log(rn^/M^), which have to be resummed 
in order to get a reasonable result. Re- 
summing all the terms of the form (a s ) p 
(a s \og(m b /M)) n (with M = m t or M = 
mw, n = 0,1, 2, ... ) for fixed p corresponds 
for p = to leading- log (LL), for p = 1 
to next-to-leading-log (NLL), and for p = 2 
to next-to-next-to-leading-log (NNLL) preci- 
sion. The previous NLL prediction, based 
on the original QCD calculations of several 

groups 14,15,16,17,18,19,20,21,22,23,2^ had an ad _ 

ditional charm mass renormalization scheme 
ambiguity, first analysed in Ref. 25 . For an 
energy cut E 1 > 1.6 GeV it reads 27 

B{B -> V s7 ) N ll = 
= (3.6lt£; 24 ± 0.02 ± 0.25 ± 0.15) x 10~ 4 (2) 

where the errors are due to the charm scheme 
dependence, CKM input, further parametric 
dependences, and to perturbative scale de- 
pendence. The dominant uncertainty related 
to the definition of m c was taken into ac- 
count by varying m c /mb in the conservative 
range 0.18 < m c /mb < 0.31, which covers 



both, the pole mass value (with its numeri- 
cal error) and the running mass value m c (/z c ) 
with /i c G [m c ,TOfc]; for the central value 
m c /rrib = 0.23 was used 27 . However, the 
renormalization scheme for m c is an NNLL 
issue. It was shown that a complete NNLL 
calculation reduces this large uncertainty at 
least by a factor of 2 28 . 

Within a global effort, such a NNLL cal- 
culation was quite recently finalized 29 . The 
calculational steps were performed by vari- 
ous groups 31,32,33,38,35,34,36,37,30,39 Qnc cru- 
cial piece is the calculation of the three-loop 
matrix elements of the four-quark operators, 
which was first made within the so-called 
large- /?o approximation 30 . A calculation 
that goes beyond this approximation by em- 
ploying an interpolation in the charm quark 
mass m c from m c > nib to the physical m c 
value has just been completed 39 . It is that 
part of the NNLL calculation where there is 
still space for improvement. 

All those results lead to the first estimate 
of the B — * X s ^f branching ratio to NNLL 
precision. It reads for a photon energy cut 
£L > 1.0 GeV 



29. 



B{B -» V s7 ) NNLL = 

(3.27 ± 0.23) x 10- 4 . (3) 

The overall uncertainty consists of non- 
perturbative (5%), parametric (3%), higher- 
order (3%) and m c -interpolation ambiguity 
(3%), which have been added in quadrature. 
For higher photon energy cut we have the 
following numerical fit: 

(4) 

where y = E /(1.0GeV) - 1. This formula 
coincides with the NNLL results up to ±0.1% 
for E Q € [1.0, 1.6] GeV. The error is prac- 
tically Efl-independent in this range. For 
E 1 > 1.6 GeV the NNLL prediction reads 29 : 



B{B -» V s7 ) NNLL = 

(3.15 ±0.23) x 10" 4 . 



(5) 
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Fig. 2. Renormalization-scale dependence of B(B — > 
X s 7) in units 10~ 4 at the LL (dotted lines), NLL 
(dashed lines) and NNLL (solid lines). The plots de- 
scribe the dependence on the matching scale /in, the 
low-energy scale fn,, and the charm mass renormaliza- 
tion scale fi c . 



Compared with the HFAG average, given in 
Eq.(l), the NNLL prediction is 1.2a below 
the experimental data (see Fig. 1 26 ). 

The reduction of the renormalization- 
scale dependence at the NNLL is clearly seen 
in Fig. 2. The most important effect occurs 
for the charm mass MS renormalization scale 
He that was the main source of uncertainty 
at the NLL. The current uncertainty of ±3% 
due to higher-order (O(afj) effects is esti- 
mated from the NNLL curves in Fig. 2. The 
reduction factor of the perturbative error is 



more than a factor 3. The central value of 
the NNLL prediction is based on the choices 
Mb = 2.5GcV and fj, c = 1.5 GeV. 

There are some perturbative NNLL cor- 
rections which arc not included yet in the 
present NNLL estimate, but are expected to 
be smaller than the current uncertainty: the 
virtual- and bremsstrahlung contributions to 
the (07, Os) and (0s, 0s) interferences at or- 
der a 2 , the NNLL bremsstrahlung contribu- 
tions in the large /^-approximation beyond 
the (07,07) interference term (which are al- 
ready available 40 ), the four- loop mixing of 
the four-quark operators into the operator 
0§ (see recent work 38 ), the exact m c depen- 
dence of the various matrix elements beyond 
the large /3n approximation (see ) and per- 
turbative logarithmically-enhanced cut ef- 
fects (see discussion below and 60 ). 

Nevertheless, the final result includes 
subdominant contributions such as the per- 
turbative electroweak two-loop corrections 
of order -3.7% 42 < 24 < 4 3,44 and thc nQn _ 

perturbative corrections scaling with 1/m 2 
or 1/to 2 of order +1% and +3% respec- 
tively 45,48,49,50,51 _ 

It is well known, that the local oper- 
ator product expansion (OPE) for the de- 
cay B — ► X s "f has certain limitations if one 
takes into account other operators than the 
leading 07, as was already shown within thc 
analysis of 1/m 2 power corrections. The ad- 
ditional error of 5% in the NNLL predic- 
tion corresponds to non-perturbative correc- 
tions, which scale with a s A/mb- Quite re- 
cently, a specific piece of the additional non- 
perturbative corrections was estimated 52 . 
Because the overall sign of the whole effect 
is still unknown, this partial estimate is not 
included in the central value of the present 
NNLL prediction 29 . 

However, there are more subtleties. 
There is an additional sensitivity to non- 
perturbative physics, due to necessary cuts 
in the photon energy spectrum to suppress 
the background from other B decays (see 



1 



Fig. 3). This leads to a breakdown of the 
local OPE, which can be cured by partial re- 
summation of these effects to all orders into 
a non-perturbative shape- function 53 > 54 > 55 . 
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Fig. 3. Cut in the photon energy spectrum 

Those shape-function effects are taken 
into account in the experimental value by 
the HFAG 13 and the corresponding theo- 
retical uncertainties due to this model de- 
pendence is reflected in the extrapolation er- 
ror in the experimental number quoted above 
in Eq. (1). Here, one should keep in mind 
that the experimental energy cuts in the 
last experiments are at 1.8 GeV or 1.9 GeV 
or even higher. The extrapolation down to 
1.6 GeV is done using three different theo- 
retical schemes to calculate the extrapolation 
factor 56,24,57,58 an( j averaging those results. 

Moreover, it was argued that a cut 
around 1.6 GeV might not guarantee that 
a theoretical description in terms of a local 
OPE is sufficient because of the sensitivity to 
the scale A = m b ~2E~ f 59 . A multiscale OPE 
with three short-distance scales m&, \fm~bA, 
and A was proposed to connect the shape 
function and the local OPE region. Quite re- 
cently, such additional cutoff-related effects 
were numerically estimated using (model- 
independent) SCET methods 60 > 61 > 62 . Those 
perturbative effects due to the additional 
scale are negligible for 1.0 GeV but lead to 
an effect of order 3% at 1.6 GeV 60 . The 
size of these effects at 1.6 GeV is at the same 



level as the 3% higher-order uncertainty in 
the present NNLL prediction. It is suggestive 
that in the future those additional perturba- 
tive cut effects get analysed and combined 
together with those already included in the 
experimental average of the HFAG. 

There are also other claims for non- 
negligible cut effects at 1.6 GeV 63 which, 
however, are based on models of the non- 
perturbative shape function. Moreover, 
there is an alternative approach to the cut 
effects in the photon energy spectrum based 
on dressed gluon exponentiation and incor- 
porating Sudakov and renormalon resumma- 
tions 64 ' 65 . It should be emphasized that 
the higher predictive power of this approach 
is related in part to the assumption that 
non-perturbative power corrections associ- 
ated with the shape function follow the pat- 
tern of ambiguities present in the perturba- 
tive calculation 66 . 



3. B -> X s £+£~ 

In comparison to the B — > X s ~{, the inclu- 
sive B — ► X s £ + £~ decay presents a comple- 
mentary and also more complex test of the 
SM. The decay B -> X s l + l~ is particularly 
attractive because of kinematic observables 
such as the invariant dilepton mass spectrum 
and the forward-backward (FB) asymmetry. 
These observables are dominated by pertur- 
bative contributions if the cc resonances that 
show up as large peaks in the dilepton in- 
variant mass spectrum are removed by ap- 
propriate kinematic cuts (see Fig. 4). In the 
'perturbative s = q 2 / m 2 -windows' , namely 
in the low-dilepton-mass region 1 GeV < 
q 2 < 6 GeV, and also in the high-dilepton- 
mass region with q 2 > 14.4 GeV, theoretical 
predictions for the invariant mass spectrum 
are dominated by the perturbative contribu- 
tions; a theoretical precision of order 10% is 
in principle possible. Regarding the choice of 
precise cuts in the dilepton mass spectrum, 
it is important that theory and experiment 
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can be compared using the same energy cuts 
and any kind of extrapolation is avoided. 




Fig. 4. Schematic dilepton-mass dependence of the 
differential branching ratio d/ds BR(B — > X 3 i+£~) in 
units 10 — 5 (the dashed line corresponds to the pertur- 
bative contribution). 

The recently calculated NNLL contri- 
butions 67,68,69,70,71,72 have s i gmncant ly im- 
proved the sensitivity of the inclusive B — > 
X s £ + t~ decay in testing extensions of the 
SM in the sector of flavour dynamics, in par- 
ticular the value of the dilcpton invariant 
mass <7q, for which the differential forward- 
backward asymmetry vanishes, is one of the 
most precise predictions in flavour physics 
with a theoretical uncertainty of order 5%. 

Also non-perturbative corrections scal- 
ing with \jm\ or l/m 2 c 47,75,51,76,77,78 are 

taken into account in the present NNLL 
predictions. The unknown non-perturbative 
corrections which scale with a s A/mb are less 
important than in the case of the decay B — » 

Recently, further refinements were pre- 
sented such as the NLL QED two-loop cor- 
rections to the Wilson coefficients of order 
2% 72 . Furthermore, it was shown that in 
the QED two-loop corrections to matrix ele- 
ments large collinear logarithms of the form 
log(mb/mi cpt on) survive integration if only a 
restricted part of the dilepton mass spectrum 
is considered. This adds another contribu- 
tion of order +2% in the low-q 2 region 73 . 

A recent update of the dilepton mass 
spectrum, integrated over the low dilepton 



invariant mass region in the muonic case, 
leads to 73 

B(B -> A s /iV~) = (1-59 ± 0.11) x 1CT 6 , 

(6) 

where the error includes the parametric and 
perturbative uncertainties only. For B(B — > 
X s e + e~), in the current BaBar and Belle se- 
tups, the logarithm of the lepton mass gets 
replaced by angular-cut parameters and the 
integrated branching ratio for the electrons is 
expected to be close to that for the muons. 
The analogous update of the other NNLL 
predictions will be presented in a forthcom- 

74 

mg paper . 

There are further subtleties, which again 
lead to larger theoretical uncertainties. In 
the high-s region, one encounters the break- 
down of the heavy-mass expansion at the 
endpoint; while the partonic contribution 
vanishes in the end-point, the 1/m 2 correc- 
tions tend towards a non-zero value. In con- 
trast to the endpoint region of the photon 
energy spectrum in the B — > X s j decay, no 
partial all-orders resummation into a shape 
function is possible here. However, for an 
integrated high-s spectrum R(s) an effective 
expansion is found in inverse powers of 

mf = m b x (1 - Vsmin) (7) 

rather than mi,. The expansion converges 
less rapidly, depending on the lower dilepton- 
mass cut s m ; n 69 . 

A hadronic invariant-mass cut is im- 
posed in the present experiments (Babar: 
m x < 1.8 GeV, Belle: m x < 2.0 GeV) 
in order to eliminate the background such 
as b — > c(— > se + v)erD = b — > se + e~ + 
missing energy. The high-dilepton mass re- 
gion is not affected by this cut and in the 
low-dilepton mass region the kinematics with 
a jet-like X s and m x < mf,AQCD implies the 
relevance of the shape function. A recent 
SCET analysis shows that using the univer- 
sality of jet and shape functions the 10 — 30% 
reduction of the dilepton mass spectrum can 
be accurately computed using the B — > X s j 
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shape function. Nevertheless effects of sub- 
leading shape functions lead to an additional 
uncertainty of 5% 79 > 80 . 

It is well-known, that the measurement 
of the dilepton mass spectrum, the zero of 
the forward-backward asymmetry, and the 
B — ► X s -y branching ratio allows to fix the 
magnitude and sign of all relevant Wilson 
coefficients within the SM. In view of the 
fact that at present only restricted data sets 
are available, it was recently proposed to fo- 
cus on quantities which are integrated (over 
q 2 ) ; besides the total rate and the integrated 
forward-bachward asymmetry it was shown 
that the third angular decomposition within 
the b — > s£ + £~ mode is sensitive to a different 
combination of Wilson coefficients 81 . 

4. Exclusive b — > s transitions 

The corresponding rare exclusive decays, 
such as B — > K*j, B — > K* 'fi + fi~ or also 
B s — ► (fry, B s — ► (f>n + [i~, are well-accessible 
at the forthcoming LHCb experiment. In 
contrast to the measurement of the branch- 
ing ratios, measurements of CP, forward- 
backward, and isospin asymmetries are less 
sensitive to hadronic uncertainties. 

For example, the value of the dilepton in- 
variant mass in B — ► K*fi + y~, q 2 , for which 
the differential forward-backward asymme- 
try vanishes, can be predicted in quite a clean 
way. In the QCD factorization approach, 
at leading order in Aqcd/^6, the value of 
<7o is free from hadronic uncertainties at or- 
der a°, a dependence on the soft form fac- 
tor £j_ and the light-cone wave functions of 
the B and K* mesons appear at order a s . 
The latter contribution, calculated within 
the QCD factorization approach, leads to a 
large shift (see 82 > 83 ' 84 ). Nevertheless, there 
is the well-known issue of power corrections 
(A-QCD/nT-b) within the QCD factorization 
approach which increases the theoretical un- 
certainty. 

An extension of the QCD factorization 



formula to the non-resonant decay B — > 
Ktt£ + £~ with an energetic Ktt pair and also 
with an energetic kaon and a soft pion was 
presented 85 . Here one relies on the fact that 
the forward-backward asymmetry is due to 
the interference of the helicity J = 1 + ,1~ 
amplitudes induced by the b — > s current. 
So it seems that no angular analysis is nec- 
essary to disentangle vector and tensor fi- 
nal states; however, the dependence of the 
non-perturbative input functions on the kine- 
matic variables might differ for the two cases. 
This suggests that a restriction to the res- 
onant states only is still the theoretically 
cleanest option. 

There are also certain transversity am- 
plitudes in B — > K*/i + (i~, in which the 
hadronic formfactors also cancel out at lead- 
ing order. Thus, such observables are rather 
insensitive to hadronic uncertainties, but 
highly sensitive to non-standard chiral struc- 
tures of the b — > s current 86 ' 87 . 

Quite recently, branching ratios, isospin 
and CP asymmetries in exclusive radiative 
decays like B — > K*j and B s — -> (f>-f were 
estimated combining QCD factorization re- 
sults with QCD sum rule estimates of power 
corrections, namely long-distance contribu- 
tions due to photon and soft-gluon emis- 
son from quark loops 88 . Particulary, this 
leads to an estimate of the time-dependent 
CP asymmetry in B° -> K* -/ of S = 
— 0.022±0.015±g 01 88 . The contribution due 
to soft gluon emssion is estimated within the 
QCD sum rule approach to be very small, 
gsgluon = _o.005±0.01, while a conservative 
dimensional estimate of this contribution due 
to a nonlocal SCET operator series leads to 
\S s 9 luon \ w 0.06 89 < 90 . Furthermore, one finds 
a larger time-dependent CP asymmetry of 
around 10% within the inclusive mode 89 ; 
however, the SCET estimate shows that the 
expansion parameter is Aqcd/Q where Q is 
the kinetic energy of the hadronic part, while 
there is no contribution at leading order. 
Therefore, the effect is expected to be larger 
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for larger invariant hadronic mass, thus, the 
K* mode has to have the smallest effect, be- 
low the 'average' 10% 89 . 

5. Comments on the so-called 
B — > Kir puzzle 

The B — > Kir modes are well known for be- 
ing sensitive to new electroweak b — > s pen- 
guins beyond the SM 91 - 92 . The data on 
CP-averaged Kir branching ratios can be ex- 
pressed in terms of three ratios: 

t b+ B[B° -> TT-K+] + B[B° -► ir+K-} 



R 



Rn. — 



t b o B[B+ 
1B[B° 



> TT+K } 
7T-K+] 



B[B° 



> TT-K } 
7T+K-] 



R r = 2 



2 B[B° 

mi - 



> ir°K } + B[B° 
ir°K+]+B\BT 



ir°K } 

► n°K-] 



B[B+ ^-K+K*]+B[B- d 



-K°] 



The actual data presented at ICHEP06 
read 94 ' 95 ' 96 



R = 0.921°;°^ = LOO^iic = 1.10^ 
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Fig. 5. Constraint in the (p, if) plane induced by the 
tttt, Kit, K K data compared with the standard CKM 
fit. 



One should emphasize that in the pre- 
vious analyses the radiative electromagnetic 



corrections to charged particles in the fi- 
nal state were not taken into account, as 
was emphasized in the past (see for exam- 
ple Ref. 104 ). These corrections, worked out 
in Ref. 105 , are now properly included in the 
analysis of both experiments. 

The present data is compatible with 
the approximate sum rule proposed in 
Rcfs. 97 . 98 ." ) which leads to the prediction 
R c rs R n . The data is also in agreement with 
the available SM approaches to these data 
based on QCD factorization techniques and 
on SU(3)f symmetry assumptions 102 > 103 . 
For example, the BBNS predictions, based 
on the QCD factorization approach 10 °. 101 ; 
are 



R 



1.16: 



0.19 



Rc 



1.15 



-0.19 
-0.17 



The latest SU(3)f analysis of the 
CKM fitter group 106 , includes all available 
7T7T, Kir, KK modes; also so-called annihi- 
lation/exchange topologies and factorizablc 
SU(3)f breaking are taken into account. 

As shown in Fig. 5 106 the constraint in 
the (p, fj) plane induced by these data im- 
plies that the compatibility with the SU(3) 
and SM hypothesis is very good (the so- 
called pValuc of that SM analysis is of order 
30 - 40%). But the xLn is not always the 
best measure of the compatibility of the data 
with the theory. Among the main contribu- 
tions to the x 2 there is the ratio R n and the 
CP asymmetry S(K° s ir°). The latter quan- 
tity is defined via 

a CP {K° s ir°}(t) = (8) 
S(K^ir°)sm(Am d t) - C(K^ir°)cos(Am d t). 

Both observables are very sensitive to new 
electroweak penguins. After removing them 
from the global fit, Fig. 6 106 shows the com- 
parison of the indirect fit (2a contour) , with 
p, fj from the CKM fit and all other avail- 
able modes, with the direct measurements 
(lcr band) using the new data. The indi- 
rect prediction for R n is now in good agree- 
ment with the direct measurement. There is 
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Fig. 6. Comparison of direct measurements of R„ and 
S(K^ir°) (la band) with the indirect fit (2cr contour). 



still a small tension in the case of the observ- 
able S(KgTT°) which is at present not really 
'puzzling' from the statistical point of view. 
Fig. 6 shows that these two quantities are not 
correlated, so that possible deviations from 
the SM values could have a completely differ- 
ent origin. These findings were just recently 
confirmed in a similar approach 107 . 

Because of the large non-factorizable 
contributions identified in the B — > 7T7t 
channel, large non-factorizable SU(3)f - 
or isospin-violating QCD and QED effects 
within the SM cannot be ruled out at 
present 108 . Future data from the B facto- 
ries and LHCb will clarify the situation com- 
pletely. There will be up to 38 measured ob- 
servables depending on the same 13+2 theo- 
retical parameters. This will allow for the 
study of SU(3) breaking and new-physics 
effects. 
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